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PROGRAM SUMMARY 

Title of the software package: ODEtools. 
Catalogue number: (supplied by Elsevier) 

Software obtainable from: CPC Program Library, Queen's University of Belfast, N. Ireland (see application form 
in this issue) 

Licensing provisions: none 

Operating systems under which the program has been tested: UNIX systems, Macintosh, DOS (AT 386, 486 and 
Pentium based) systems, DEC VMS, IBM CMS. 

Programming language used: Maple V Release 3/4. 

Memory required to execute with typical data: 16 Megabytes. 

No. of lines in distributed program, excluding On-Line Help, etc.: 10159. 

Keywords: 1 st /2 nd order ordinary differential equations, symmetry methods, symbolic computation. 
Nature of mathematical problem 

Analytical solving of 1 st and 2 nd order ordinary differential equations using symmetry methods, and the inverse 
problem; that is: given a set of point and/or dynamical symmetries, to find the most general invariant 1 st or 
2 nd order ODE. 

Methods of solution 

Computer algebra implementation of Lie group symmetry methods. 
Restrictions concerning the complexity of the problem 

Besides the inherent restrictions of the method (there is as yet no general scheme for solving the associated PDE 
for the coefficients of the infinitesimal symmetry generator), the present implementation does not work with 
systems of ODEs nor with ODEs of differential order higher than two. 

Typical running time 

This depends strongly on the ODE to be solved. For the case of first order ODEs, it usually takes from a few 
seconds to 1 or 2 minutes. In the tests we ran with the first 500 1 st order ODEs from Kamke's book Jl]], the 
average times were: 8 sec. for a solved ODE and 15 sec. for an unsolved one, using a Pentium 200 with 64 
Mb. RAM, on a Windows 95 platform. In the case of second order ODEs, the average times for the non-linear 
2 nd order examples of Kamke's Book were 35 seconds for a solved ODE and 50 seconds for an unsolved one. The 
tests were run using the Maple version under development, but almost equivalent results are obtained using the 
available Maple R4 and R3 (the code presented in this work runs in all these versions). 

Unusual features of the program 

The ODE-solver here presented is an implementation of all the steps of the symmetry method solving scheme; 
that is, the command receives an ODE, and when successful it directly returns a closed form solution for the 
undetermined function. Also, this solver permits the user to optionally participate in the solving process by giving 
advice concerning the functional form for the coefficients of the infinitesimal symmetry generator (infinitesimals) . 
Many of the intermediate steps of the symmetry scheme are available as user-level commands too. Using the 
package's commands, it is then possible to obtain the infinitesimals, the related canonical coordinates, the finite 
form of the related group transformation equations, etc. Routines for testing the returned results, especially 
when they come in implicit form, are also provided. Special efforts were put in commands for solving the inverse 
problem too; that is, commands returning the most general 1 st or 2 nd order ODE simultaneously invariant under 
given symmetries. One of the striking new features of the package -related to 2 nd order ODEs- is its ability to deal 
with dynamical symmetries, both in finding them and in using them in the integration procedures. Finally, the 
package also includes a command for classifying ODEs, optionally popping up Help pages based on Kamke's advice 
for solving them, facilitating the study of a given ODE and the use of the package with pedagogical purposes. 
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LONG WRITE-UP 



Introduction 

In a previous work , we presented an implementation of Lie symmetry methods (SM) for solving first 
order ODEs. The key idea of that work was to prepare routines for finding particular solutions for the 
PDE determining the coefficients of the infinitesimal symmetry generators (infinitesimals), as well as 
providing extra routines for the user's input of functional form ansatze when the default routines fail. 
This approach is presently solving 85% of Kamke's examples using only the defaults, apart from being a 
concrete way to tackle non- classifiable 1 st order ODEs, for which the standard computer algebra solvers 
usually fail. 

The same idea can be implemented for 2 nd order ODEs too. To understand the motivation, we recall 
that most implementations of SM for high order ODEs are based on the setup and solving of the so-called 
determining equations for the infinitesimals -an overdetermined system of PDEs-, which arise when we 
assume we are interested only in point symmetries ||. However, this poses a limitation on the 2 nd order 
ODEs that can be dealt with, since only a restricted subset of them have point symmetries. On the other 
hand, all ODEs have infinite dynamical symmetries |Q, |5|, which arise as particular solutions for a single 
linear 2 nd order PDE for the infinitesimals. 

This paper then presents the implementation of the ideas of our previous work to tackle 2 nd order 
ODEs. One of the tricky things related to the use of SM to solve n th order ODEs is that the knowledge of 
n symmetries does not directly reduce the problem to a line integral as in the 1 st order case. Moreover, the 
alternatives found in the literature for constructing the solution departing from dynamical symmetries are 
few. We then extended some of the standard integration methods for point symmetries, and implemented 
them as routines for dynamical symmetries too. 

As a second issue, we invested in the research design of the package, extending both the number and 
the capabilities of extra user-lever routines related to the intermediate steps of the symmetry scheme. 
Worth mentioning are a routine for finding the most general 1 st /2 nd order ODE simultaneously invariant 
under many point or dynamical symmetries, and a routine returning the symmetries of an unknown ODE, 
given its solution. 

As a third issue, we invested in augmenting the pedagogical potential of the package by extending the 
classification capabilities of the odeadvisor command to work with most of the standard classifications 
for 2 nd order ODEs (see §). 

The exposition is organized as follows. In sec. |l|, the SM scheme for solving 2 nd order ODEs is briefly 
reviewed. In sec. ||, a compact table-summary of all ODEtools routines and a detailed description of 
the ODE-solver are presented. Sec. |^ briefly illustrates the extension of the methods presented in Q 
for finding the infinitesimals. In sec. |], we comment about the methods implemented for integrating 
ODEs from the knowledge of their symmetries, focusing on the case in which these symmetries are of 
dynamical type. Sec. |^ displays the results of testing the package with the non-linear 2 nd order examples 
of Kamke's book as well as an update of the results obtained for the Kamke's first 500 1 st order ODEs. 
In sec. | the main differences between ODEtools and other existing packages for symmetry analysis of 
ODEs are highlighted. Finally, the conclusions contain a brief discussion about this work and its possible 
extensions. 



1 Symmetry methods for 2 nd order ODEs 

Generally speaking, the key point of Lie's method for solving ODEs is that the knowledge of a (Lie) group 
of transformations which leaves a given ODE invariant may help in reducing the problem of finding its 
solution to quadratures J?], |^, [)) . Aside from the subtleties which arise when considering different cases, 
we can summarize the computational task of using SM for solving a given 2 nd order ODE, say, 



-7-2 = ®( x > y> -r< 

ax dx 

as the finding of the infinitesimals of a one-parameter Lie group which leaves Eq.(Q) invariant, i.e., a pair 
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of functions^ {f(y, x, yx), T](y, x, yx)} satisfying 
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followed by the integration of the ODE by either: 

a) determining the differential invariants of order and 1 of the symmetry group; 

b) determining first integrals ip(x,y,yi) = ipo (from the knowledge of £ and 77); 

c) determining the canonical coordinates, say {r, s(r)}, of the associated Lie group. 

A first look at the symmetry scheme may lead to the conclusion that finding solutions to Eq.(|^) would 
be much more difficult than solving the original Eq.(Q). However, what we are really looking for is a 
particular solution to Eq. (Q) , and in many cases this particular solution is the only thing one can actually 
obtain. As an example, consider 

This 2 nd order ODE is non-linear and it doesn't match any pattern for which we know the solving method a 
priori; standard classification based ODE-solvers then fail when trying to solve it. However, a polynomial 
ansatz for the infinitesimals (here made by symgen, the routine for determining the symmetries) rapidly 
leads to 3 particular solutions to Eq.(^)^: 

> symgen ( ode ) ; # input = DDE, output = pairs of infinitesimals 



[£ = 0,ry = x], [£ = x,ri = y], [( = x 2 ,r] = yx] (4) 

Passing the ODE directly to odsolve (the solver) , it will internally call symgen and use the result above 
to solve the ODE as follows: 

> odsolve (ode) ; 

y = {\n{x) - ln(l + xC x ) + C 2 ) x (5) 

What is amazing, and characteristic of symmetry methods is that if we change (ji^x — y^j /xtoF (f^x — y 
in Eq.(^|), where F is an arbitrary function of its argument, the first two symmetries of Eq.(Q) will remain 
valid and the solving scheme will succeed as well: 

dx 2 x \\dx} x} ^ ^ 

> odsolve (ode) ; 
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In what follows, y\ = -X. 
In what follows, the input can be recognized by the Maple prompt >. 
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y = 



RootOf ln(a:) + d 



1 , \dx 

(_o-F(_o)) j x 



C 2 



(7) 



The answer above is expressed using RootOf, and the inner integral uses the new intat command^. 
This kind of general answer can be interpreted as a mapping in that, given F, it returns the answer after 
calculating the roots of the resulting expression; for example, the ODE Eq.@ is a particular case of 
Eq.(^|), and its answer Eq.(^|) is what one would obtain introducing F : u — > u 2 in Eq.(^). 

Also, due to the fact that we are just looking for particular solutions to Eq.(@), symmetry methods 
can be an alternative in solving linear ODEs too. Consider, for instance: 



dx 2 i i dx J x 



(8) 



This ODE does not match a known pattern related to special functions, nor can it be solved using the 
standard schemes for rational or exponential solutions, due to the presence of an arbitrary function F{x). 
However, its invariance under [£ = 0, r\ = x] and [£ = 0, r\ = y] can be easily determined by a polynomial 
ansatz, from which its solution follows straightforwardly: 



> odsolve (ode) ; 



y = x 



xF(x) 



-d.r \ \ 

dx + C 2 

J 



(9) 



2 The package's commands 

A detailed description of the package's commands, with examples and explanations concerning their 
calling sequences, is found in the On-Line Help, and is already present in [Q. Therefore, we have restricted 
this section to a brief table-summary and a detailed description only of the solver, odsolve. Some 
input/output examples can be seen in sec. |^ and ||. 



2 . 1 Summary 

A compact summary of the commands of the package is as follows: 



Command 


Purpose: 


odsolve 


solves ODEs using the symmetry method scheme 


intfactor 


looks for an integrating factor for first order ODEs 


canoni 


looks for a pair of canonical coordinates of a given Lie group 


eta_k 


returns the k-extended infinitesimal 


transinv 


looks for the finite group transformation 


infgen 


returns the k-extended symmetry generator as an operator 


symgen 


looks for pairs of infinitesimals 


equinv 


looks for the most general ODE invariant under a given set of symmetries 


buildsym 


looks for the infinitesimals given the solution of an ODE 


odepde 


returns the PDE for the infinitesimals 


odetest 


tests explicit/implicit results obtained by ODE-solvers 


symtest 


tests a given symmetry w.r.t a given ODE 


odeadvisor 


classifies 1 st /2 nd order ODEs and pop up related Help-pages 



Table 1. Summary of the ODEtools commands 



4 intat is a command of the last version of PDEtools [[UJ , and represents an integral evaluated at a point -analogous to 
a derivative evaluated at a point, intat displays the evaluation point as an upper limit of integration. 
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2.2 Description 

Command name: odsolve 

Feature: 1 st and 2 nd order ODE-solver based on symmetry methods 
Calling sequence: 

> odsolve (ode) ; 

> odsolve(ode, y(x), way=xxx, HINT= [ [exprl , expr2],..], int_scheme) ; 

Parameters: 

ode - a 1 st or 2 nd order ODE. 

y(x) - the dependent variable (required when not obvious). 

way=xxx - optional, forces the use of only one (xxx) of the 8 internal 

algorithms {abacol, 2, 3, 4, 5, 6, abaco2, pdsolve} for deter- 
mining the coefficients of the infinitesimal symmetry generator 
(infinitesimals). 

HINT = [el,e2] - optional, el and e2, indicate possible functional forms for the 

infinitesimals. 

HINT= [ [el , e2] , [e3,e4] , . .] - optional, a list of hints for the infinitesimals. 

int_scheme - optional, one of: fat, can, can2, gon, gon2, dif. 

Optional parameters can be given alone or in conjunction, and in any order. 

Synopsis: 

Given a 1 st or 2 nd order ODE, odsolve's main goal is to solve it in two steps: first, determine pairs 
of infinitesimals of 1-parameter symmetry groups which leave the ODE invariant, and then use these 
infinitesimals to integrate the ODE. 

To determine the infinitesimals, odsolve makes calls to symgen, another command of the package. 
To integrate the ODE using these infinitesimals, odsolve has seven schemes, almost all of them explained 
in connection with point symmetries in |?], |§| : 

1. building an integrating factor (fat, only for 1st order ODEs) 

2. reducing the ODE to a quadrature using the canonical coordinates of that group (can) 

3. reducing a 2 nd order ODE to a quadrature at once, using 2 pairs of infinitesimals forming a 2-D 
subalgebra (can2) 

4. reducing a 2 nd order ODE to a quadrature, constructively, using a normal form of the generator in 
the space of first integrals (gon) 

5. reducing a 2 nd order ODE to a quadrature at once, using 2 pairs of infinitesimals (a 2-D subalgebra) 
and normal forms of generators in the space of first integrals (gon2) 

6. using differential invariants constructively (dif). 

7. solving a 2nd order ODE at once, using 3 pairs of infinitesimals, when no two of them can be used 
to form a 2-D subalgebra. 

The integration schemes ([j]) and (Q) are used with 1 st order ODEs, while schemes (|^) to (^) work 
with 2 nd order ODEs. The integration schemes (f|), (||) and (||) work with dynamical symmetries too. 
odsolve does not classify the ODE before tackling it and is mainly concerned with non-classifiable ODEs 
for which the standard Maple dsolve fails. By default, odsolve starts off trying to isolate the derivative 
in the given ODE, then sequentially uses subsets of the algorithms of symgen to try to determine the 
infinitesimals, and finally sequentially tries the integration schemes mentioned above. The default order 
for trying these schemes is 

• 1 st order ODEs: can, fat 

• 2 nd order ODEs: gon, can2, gon2, can, dif 

When odsolve succeeds in solving the ODE, it returns, in order of preference: 
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• an explicit closed form solution; 

• an implicit closed form solution; 

• a parametric solution (default strategy for 1 st order ODEs when the derivative cannot be isolated). 
All these defaults can be changed by the user; the main options she/he has are: 

• To request the use of only one of the algorithms for determining the infinitesimals (way=xxx option; 
different algorithms may lead to different symmetries for one and the same ODE, sometimes making 
the integration step easier). 

• To enforce the use, in a specific order, of only one or more of the alternative schemes for integrating 
a given ODE after finding the infinitesimals (fat, can, can2, gon, gon2, and dif optional arguments, 
useful to select the best integration strategy for each case). 

• To indicate a possible functional form for the infinitesimals (HINT=xxx option). This option is 
valuable when the solver fails, or to study the possible connection between the algebraic pattern of 
the given ODE and that of the symmetry generators. 

A brief description of how the HINT=xxx option can be used is as follows: 

• HINT= [el , e2] , indicates to the solver that it should take el and e2 as the infinitesimals and 
determine the form of (a maximum of two) indeterminate functions possibly contained in el and/or 
e2, such as to solve the problem. 

• HINT= [ [el , e2] , [e3 , e4] ,...], where [el , e2] is any of the above. 

• HINT=parametric, indicates to the solver that it should only look for a parametric solution for the 
given (1 st order) ODE. 

Finally, there are three global variables managing the solving process, which are automatically set by 
internal routines but can also be assigned by the user, as desired. They are dgun, ngun, sgun, for setting, 
respectively, the maximum degree of polynomials entering some of the ansatze for the infinitesimals, the 
maximum number of subproblems into which the original ODE should be mapped, and the maximum 
size permitted for such subproblems. The dgun variable is automatically set each time symgen is called, 
according to the given ODE, whereas, by default, ngun and sgun have their values assigned to 1. Increasing 
the value of dgun usually helps, especially in the case of polynomial ODEs; but, although the user- 
assigning of the ngun or sgun variables might increase the efficacy of the algorithms, each increase of one 
unit can slow down the solving process geometrically. 



3 Finding the infinitesimals 

In the context of SM, the infinitesimals we are looking for are solutions of Eq.(Q), which is linear in the 
functions (£, ?y) and their derivatives. The command which looks for these pairs of infinitesimals for a 
given ODE is symgen, and the version here presented is an extension to 2 nd order ODEs of the schemes 
presented in 

3.1 The symgen subroutines 

The key idea underlying this extension is that to look for dynamical symmetries for 2 nd order ODEs 
is mainly equivalent to looking for point symmetries for 1 st order ODEs; except that the routines are 
now going to look for solutions involving (x, y, y±) instead of only (x,y). More specifically, we adapted 
the previous symgen/... subroutines, whose main purpose is to look for particular solutions to linear 
problems, in order to look for such particular solutions considering y\ as a new variable in equal footing 
as x and y. An explanation with details about how these symgen/... routines work can be found in Q; 
so, we here restricted the discussion to some examples illustrating the type of results which can now be 
obtained. 
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The I s * algorithm, abacoi, typically looks for infinitesimals of the form [£ — Q,j] — F(q)] and [£ 
F(q), 77 = 0] , where q is one of (x, y, y\). 
Example: Kamke's 2 nd order non-linear ODE 206 



(a 2 - x 2 ) (a 2 - t, 2 ) g + (a 2 - x 2 ) „ (g) - * (a 2 - y 
> symgen (ode ,way=abacol) ; 



..«!,:= (<r- .a(--/r)T4-(---)//(^) 2 -x(a 2 - y 2 )g=0 (10) 



K = ' = V-^ + y 2 (l + Info + V-a'+y 2 ))], 

[i = Vx 2 - a 2 (l + hi(x + y/x 2 - a 2 )) , 77 = 0] (11) 

It is also not difficult to find the patterns of 2 nd order ODEs having the type of symmetries which abaco\ 
is prepared in principle to look for. These ODE-patterns can be obtained using the equin v ro utine, 
programmed to solve the inverse problem (mpitt=symmetries, output— invariant ODE; see sec. 3^). For 
example, in the first of the six cases mentioned above, 

> equinv([F(x) ,0] ,y(x) ,2) ; 

is the most general 2 nd order ODE invariant under [.F(a;),0], F and _F1 being arbitrary functions of its 
arguments, and Eq.(|lC|) is a particular case of this ODE family. 

The 2 nd and 3 rd algorithms look for polynomial in (x,y,yi) solutions to Eq.(|2j). 
Example: Kamke's 2 nd order non-linear ODE 181 



ode:= x 2 (x + y) - (x^ - y) = <> (1:-!) 



> symgen ( ode , way=3 ) ; 

[£ = x,r) = y], [£ = -x,ri = x], [£ = x 2 , 77 = xy] (14) 

> odsolve (ode) ; 

( - ') 

y = xe \ ' — x (15) 

The 4 th algorithm looks for rational in (x,y,yi) solutions to Eq.(||). 
Example 



ode := -—j = 2 — + 1 — (16) 

ax z x + y \ ax J ax 



> symgen ( ode ,way=4) ; 



rre 1 m re V nl u x(x + 2y) 

x + y x + y x + y 



K = -1,77 = 1], [Z = x,v = y]} (17) 



> odsolve (ode ,way=4, can2) ; 



(i 



x 2 + 2C 1 

y = ~2^T2^ (18) 

The 5 th algorithm uses a polynomial of degree two constructed from a basis of functions and algebraic 
objects, together with their derivatives, taken from the given ODE (see Q). 
Example: Kamke's 2 nd order non-linear ODE 238 

ode:= J( , + ^(£y_,( l + 4*)£ +a (, + *}*_ „_o (19) 



V dx 2 J \ dx J dx 2 \ dx J dx 



> symgen(ode ,way=5) 



4 yi + 2 x + x 3 - xy/x 4 - 8x 3 yi - IGxyi - 16 yi 2 + 16jy + \Qx 2 y 

[£ =0,11= -7— A 2 ( 20 ) 

In the example above, the square root, which appeared after isolating — ^, was viewed as an algebraic 

dx 2 

object, and the symmetry found is of dynamical type. 

The next algorithm, abaco2, looks for infinitesimals of the form [£ = F(qi),rj = Gfa)], where q\ and 
(72 are members of {x,y,yi}. 
Example 

d 2 y ( dy\ 2 e v ( 1 . x ( dy 



ode: = i? = \t) -^ i? U +lnW '^ U»» (21) 



> symgen(ode ,way=abaco2) ; 



K = ^ >7 = e y ] 



Apart from setting up the symgen/... routines to work with 2 nd order ODEs, two new algorithms 
for 1 st order ODEs were implemented. The first one is specific for Riccati-type ODEs and the second one 
is based on a direct attempt to solve the problem using the pdsolve routine from the PDEtools package. 
The routine for Riccati ODEs makes three sequential attempts to solve Eq.(H|), the key idea being to map 
the problem of finding the infinitesimals for a given Riccati ODE into that of finding particular solutions 
to a related linear 2 nd order PDE for rj(x, y). These three attempts consist of tackling this PDE by: using 
the 5 th algorithm explained above, trying separation of variables, and using the abaco\ scheme. 
Example: Kamke's 1 st order ODE 16 



> symgen(ode ,way=6) 



ode:= ^+y 2 + (xy-l)f(x) = (22) 



> odsolve (ode ,way=6) 



x 



f(x)xdx — 2 ln(a;) 

( r-^ ) 



v = - + e -^„ ; \ (24) 



dx — C\ 
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3.2 The inverse problems 



For many reasons, it appeared interesting to develop routines to solve the inverse problems too; that is: 
find the most general 1 st /2 nd order ODE simultaneously invariant under a given set of symmetries, or 
given the answer of an unknown ODE, find its symmetries (and then the corresponding invariant ODE). 
Typical applications for these routines are: we look for a description of a problem for which we only know 
the symmetries or some particular solutions; or we look for solving methods for ODEs having a given 
pattern, for which we know the solution in some cases, and we want to find the solution for the general 
case. 

The ODEtools routines related to solving the inverse problems are equinv and buildsym. The former 
is prepared to work, in principle, with an "arbitrary number" of symmetries. As a first example, consider 
the determination of the most general 2 nd order ODEs simultaneously invariant under 

X:= [[e=l,»/=l],K = x,» 7 = y],[e = x 2 , f / = y 2 ]] (25) 

> equinv (X, y(x), 2); 



dx 2 




(26) 



As an example involving dynamical symmetries, consider the determination of the 2 nd order ODE 
simultaneously invariant under 



X:= m = yi,r 1 = 0], [£ = 0,v=-}Y, (27) 

yi 



> equinv(X, y(x) , 2); 



'dy 

l y \dx 



dx \ dx I \dx 



(28) 



where -Fl is an arbitrary function of its argument. Proceeding as in the examples above, it is possible 
to determine the general ODE-patterns which can be solved by the abacoi scheme and some of those which 
can be solved by abaco2, as well as to build matching-pattern routines to match families of ODEs for 
which we know the form of the symmetries a priori. Of course, since the number of possible ODE families 
is infinite, it is not realistic to think of producing matching pattern routines for all the cases. On the 
other hand, for a variety of concrete problems which arise in different applied areas, it may be convenient 
to build these routines and have the related solving schemes at hand. Remarkably, most of the standard 
methods actually work in this manner (i.e., by matching a pattern and then applying the corresponding 
solving scheme), and these patterns actually arose from applied problems. As a benchmark for equinv, 
the routine succeeded in finding the most general invariant ODEs for all the twenty examples found in 
sec. 8.3 or Ibragimov's book on Lie methods. 

Concerning the computational scheme used in equinv, the idea is, roughly speaking, to introduce the 
given symmetries into Eq.(||), and solve it for the right-hand-side of the ODE, which is now viewed as 
"the unknown of the problem" . One of the interesting advantages of this simple approach is that we can 
consider the invariance of an ODE under point and dynamical symmetries in equal footing, and forget 
about the technicalities involved in defining differential invariants when derivatives are present in the 
transformation equations^]. 

When many symmetries are given at once, from the computational point of view, the problem is a bit 
more complicate. A brief summary of the scheme implemented is as follows: 

1. introduce the first pair of infinitesimals into Eq.(0) setting up a PDE for the right-hand-side (RHS) 
- here denominated y, y±) - of the most general explicit invariant second order ODE; 



5 This scheme works just as well to determine high order ODEs (differential order > 2) invariant under dynamical 
symmetries. 
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2. solve this PDE for <E>; the solution will involve an arbitrary function of two differential invariants, 
say F(ip 1 (x,y,y 1 ),tp 2 (x,y,y 1 )); 

3. update the symbolic value of $ introducing the result of the previous step (i.e., introducing F). 
This updated $(x,y,yi) is an explicit expression of {x, y, y{\ hereafter called (f>(x,y,yi). When 
only one symmetry is given, 4>{x, y, y\) is already the RHS of the invariant ODE we are looking for; 

4. take Eq.(^J) again and set up a new PDE using the second given pair of infinitesimals, and substitute 
$ found in PDE by its updated value <fi; 

5. we now need to solve this new PDE for the arbitrary function F introduced in step (||), but 
this is not directly possible since the unknown F depends on algebraic expressions - il>\{x, y, y\) and 
ipiix, y, y\); so introduce these differential invariants as new variables to obtain a PDE for F(tpi, ^2) 
where (ipi,ip2) are now symbol variables; 

6. solve the PDE of the previous step for F and use the inverse transformation to restore the original 
variables {x, y, y{\ in the answer; 

7. update the symbolic value of <E> with respect to the result of the previous step. The resulting ODE 
y = 4>{x, y, y\) is now simultaneously invariant under the first two pairs of infinitesimals; 

8. if more pairs of infinitesimals are given, repeat this process from step (^); 

To solve all the intermediate PDEs we use the pdsolve command of the PDEtools package |l0|, with 
which ODEtools is completely integrated^. What requires a bit more of care is the determination and 
handling of all the changes of variables involved at each step. 

Concerning the buildsym command, the idea is to build the symmetries of an unknown ODE de- 
parting from its solution. When the command succeeds it is then possible to use the equinv command 
to build the most general ODE having the symmetries of the given ODE solution. As a simple example, 
consider again Eq.(||). Passing the answer of this ODE (Eq. ([)])) to buildsym it returns two symmetries: 

> X := buildsym(ans , y(x)); 

/J-Fl (x)dx 
~ ~2 dx \ ( 29 ) 

In this case, since the ODE is linear and the symmetries are of the form [0, F{x)] , the returned symmetries 
are actually two particular solutions to Eq.(0). A less simpler example is given by: 



d 2 y x 2 dy 1 f dy^ 2 



dx 2 Jy dx 2y \dx 
An implicit solution to this ODE is given by: 



(30) 



-y 3 / 2 + ^x 3 + 2x 2 + ix-2C 1 e x -C 2 = (31) 
This ODE solution leads to the following pair of underlying symmetries 
> X := buildsym(ans , y(x)); 

X:=K = 0, r} = ^=], K = 0, r? = 4=] ( 32 ) 

In turn this information leads to the general ODE family having a solution with the same two symmetries 
of Eq.@: 



> equinv (X, y(x)); 



6 In turn, the solving of intermediate ODEs built by pdsolve to solve PDEs is carried up by ODEtools commands. 
PDEtools and ODEtools are both compiled in a single binary maple library. 
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fy 

dx 2 



-Fl (x) 

~7v~ 



dx 



1 

2j/ 



(33) 



Concerning the computational scheme implemented in buildsym, the idea is, roughly speaking, to 
use the ODE-solution to build two first integrals, from which two pairs of infinitesimals can be obtained. 
To build these pairs of infinitesimals, we note that since the first integrals, say ip and ip, are solutions of 
the linear operator 



=/ 



dx dy 



®(x,y,yi) 



df_ 

dyi 



(34) 



and from the definition of a symmetry X such that it satisfies [X, A] = XA (sec [frl) it is always possible 
to formulate the problem as: to determine two infinitesimal generators X\ and X 2 of the form [0, rj\ 
satisfying 



Xitpx 



dip 
m dy 



(l) dip 
Vi -r- = 0, 
dyi 



dip (i) dip 
dy dyi 



dp (i) dip dip (1) dip 

x 2 ipi = m~r + m ~r~ = !> x 2 p> 2 = m-r + r i2 S— = °- 

dy dyi dy dyi 

This system of equations can be solved for {rji 772}, arriving at 



Vi = 






dp 
Wl 




1 


dip 

dyi 


1 


dip 

dyi 







dip 

dyi 



dip 


dip 


- m = - 


d(p 


dtp 


dy 






dy 




dip 






dip 




dy 


dyi 




dy 


dyi 



(35) 



The determination of X\ and X 2 then amounts to the determination of ip and <p. These first integrals 
can be determined from the given solution as follows: 

1. solve the given ODE solution with respect to one of the integration constants, say C 2 , obtaining 

C 2 = F(x,y(x),C 1 ); 

2. differentiate the result of the previous item w.r.t x and solve it w.r.t Ci obtaining the first integral 
(a solution of Eq.(|34|)) C\ — p(x,y,yi); 

3. repeat the two previous items interchanging Ci <-* C 2 to obtain the second first integral C 2 = 



4 Integrating the ODE after finding the infinitesimals 

The methods implemented for integrating the ODE once the infinitesimals were found (methods (a), (b) 
and (c) mentioned in sec. are extensions to those found in 0, ^J, The extension consists of making it 
possible to use them to integrate the given ODE even when the symmetries found are of dynamical type. 
Since their use in the framework of point symmetries is rather standard, we restricted the discussion here 
to the extension of these schemes to work with dynamical symmetries. 

Concerning "why" so many integration strategies were implemented, we would like to recall that the 
knowledge of 2 symmetries of a 2 nd order ODE does not guarantee a successful integration of it since the 
different methods available do not apply in all the cases. Also, two methods which may be "equivalent" 
from the mathematical point of view, may not be equivalent for the computer algebra system, since the 
intermediate tasks involved in each integration process are different. We then preferred to implement 
all these methods, provide a default for their use, and make them available to the user by giving extra 
arguments to odsolve (see sec. ^.2[). 
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4.1 Using differential invariants 

As is well known, once a pair of infinitesimals for a given 2 nd order ODE is found, it is possible to integrate 
the ODE if one succeeds in using these infinitesimals to find the differential invariants of order and 
1 and in solving an intermediate 1 st order ODE (see for instance section 3.3.2 of H). However, if the 
symmetry is of dynamical type, the invariant of order will be a function of y x , the invariant of order 1 
will be a function of 1/2, and so on. Then, the whole scheme breaks down since, to reduce the order, the 
scheme requires that the differential invariants of order n do not depend on derivatives of order higher 
than n. 

The alternative we implemented is based on the observation that, if we transform the symmetry to 
the evolutionary forrn^ |S|, then x will always be an invariant (of order 0), so it is possible to work as 
follows: 

1. put the symmetry in the evolutionary form; 

2. use u — x as the invariant of order and find the invariant of order 1, v = v(x, y,yi)', 

3. construct the invariant of order 2: (^); 

4. change the variables and write the ODE in terms of u, v and ^ obtaining the expected reduction 
of order; 

5. solve this 1 st order ODE, and reintroduce the original variables (x, y, yi), obtaining another 1 st order 
ODE, which can be integrated straightforwardly since it has the same symmetry of the original 
2 nd order ODE. 

A tricky point in this last step is that, when the symmetry is of dynamical type (depends on to 
obtain the symmetry which holds for the 1 st order ODE one should replace all occurrences of y\ in the 
infinitesimals of the original ODE by the right-hand-side of this 1 st order ODE. 



4.2 Using first integrals 

Another way of integrating a 2 nd order ODE is to use 2 pairs of infinitesimals to build 2 first integrals, 
from which one can expect to obtain a closed form solution (see section 7.3 of @). One of the advantages 
of this method is that one always work with the original variables avoiding potential problems related 
to the introduction of changes of variables (e.g., the system fails in finding the inverse transformation 
etc.). Taking into account the way we use "dynamical symmetries" to solve the 1 st order ODE obtained 



in the reduction process (see sec. 4.1), the extension of this scheme to work with such symmetries was 
straightforward. 

Now, it is also possible to tackle a 2 nd order ODE from the knowledge of only 1 pair of infinitesimals 
(either of point or dynamical type), by looking for first integrals, exploiting similar ideas as those of the 
previous subsection; that is: 

(i) Put the known symmetry in the evolutionary form. ; 

One can then prove that there exists a first integral f(x, y, y{) which satisfies: 

A{p) = 

X(tp) = (36) 
where A is given by Eq. (|3~i|) , and X is the infinitesimal symmetry generator 

X:=f^> £(x, y, yi) |^ + 77(21, y, y x ) ^j- + r)i(x, y, y x ) 

(ii) Solve the characteristic strip of X(ip) = to obtain the invariants of order and 1: cq(x, y, y\) and 
ci(x,y,yi). 



7 Since [£ = 1, 77 = y{\ is always a (trivial) symmetry, a general symmetry can always be "gauged", that is, rewritten in 
many different forms. The evolutionary form is given by [£ = 0,r) = <3>(:r, y, 1/1)]. 
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Since the symmetry is in evolutionary form, it always happens that cq = x, and the characteristic strip is 
reduced to a single 1 st order ODE. Now, observing that the solution to X(ip) = can always be written 
as <p(co,ci), 

(iii) change the variables in A(ip) = 0, from (x,y,yi) to (co,y, Ci). 
One can then prove that A(tp) = becomes 

£ + = 

that is, a PDE in only two variables; solve it, obtaining (p{ca,c\) now simultaneously satisfying 
equations (36). Finally, 

(iv) reintroduce the original variables to obtain tp(x, y,yi); that is, a 1 st order ODE which can be 
integrated straightforwardly since it already satisfies X(ip) = 0. 



4.3 Using canonical coordinates 

It is also possible to use a pair of infinitesimals to constructively reduce the order of a 2 nd order ODE 
by introducing the canonical variables of the related invariance group. Also, if one knows two or more 
symmetries at once, one can use this information to integrate the ODE in almost only one step (see sec. 

7.4 of Q]). These two methods were implemented too, although their extension to work with dynamical 
symmetries doesn't appear to us as advantageous if compared to the extensions outlined in the previous 
subsections. 



4.4 Examples 

We present here some examples to illustrate how these integration schemes work and the different types 
of answers they may return. As the first example, consider Kamke's non- linear 2 nd order ODE 99: 

ode:= * 4 Sf + f^- y y=0 (37) 



dx 2 \ dx 

This ODE has simple polynomial point symmetries: 
> symgen(ode) ; 

[£ = 0,77 = 2], [£ = x,ri = y] (38) 



Its integration is rather fast using the method explained in sec. 4.2 
> odsolve(ode,gon) ; 

y =(-arctan(- 7 ==i===) + C7 2 )x (39) 

Remarkably, if in this case we use the integration scheme which uses the 2 symmetries to build 2 first 
integrals, the answer will not be so compact due to the manipulations required to transform the two first 
integrals tpi(x,y,yi), tf2(x,y,yi) into the solution (ps(x,y). The opposite happens with Kamke's ODE 
227: 

for which symgen finds two point symmetries, [£ = x, r] = 0] and [£ = 0, rj — y], and the more compact 
answer is returned, using special functions, by this scheme which builds two first integrals: 



> odsolve (ode ,gon2) ; 
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,.i (it^f-^Uiy^h^hU (4i) 



2 V V e 

As the next example, consider Kamke's ODE 122: 

This example involves two arbitrary functions of x and the package defaults only succeed in finding 1 
symmetry: 

> symgen(ode) ; 

[Z = 0,ri = y] (43) 

In this situation only the methods prepared to work with just 1 symmetry will work; in this example, all 
of them lead to the same result: 



ffMdx r (-ff(x)dx} Jf(x)dx N 

e g(x)e dx + e C\dx + C2 



> odsolve(ode ,gon) ; 



y = e ' (44) 

The following example, Kamke's ODE number 20, is rather critical in the sense that all the integration 
schemes fail in solving the problem, except the least elaborated one (constructive reduction of order 
introducing canonical coordinates, using only 1 symmetry at a time). Due to the presence of an arbitrary 
function of x and y(x), the answer appears in implicit form: 

S-^'K*) (45) 



> odsolve (ode , can) ; 



y 

1 - 



ln(s)-/V* d_b-C 2 = (46) 

J.b 2 + 8 F(_b)d.b + d 

As the last example, consider the case in which only 1 symmetry has been found, and the integration 
schemes only succeed in reducing the order but not in solving the resulting 1 st order ODE; e.g., Kamke's 
ODE number 225: 

d 2 y . / dy\ 2 2 ^ / 1 ( dy 



odc „ F („_! _ D(F)( „) ^ j _ {F(y) y t ^ ,_ j _ (47) 

A symmetry for this ODE is found by the fifth algorithm, 

> symgen(ode,way=5) ; 

[0,F(V)] (48) 
In these cases, an answer with the structure of the solution is returned: 

> odsolve (ode, way=5) ; 



y = RootOf Ci 



J. si (_r)eL 



F(_a) 



d-d 



&where 



J7r- Sl{ - r) = *(- r > J ' (-'))} ' [~ r = *> ~ sl ^ = F^y) fx 



(49) 
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In the structure above, we see the answer for y up to a change of variables. More specifically, this 
change of variables was not performed because the system doesn't know how to solve the displayed 
intermediate 1 st order ODEs. Although in this particular case it is impossible to go ahead -the unsolved 
1 st order ODE is the most general we can imagine - in many cases the user may be able to solve the 
reduced first order ODE by other means. Once _sl (_r) is obtained, a change of variables using the 
displayed transformation, will lead to the desired answer for y. 

5 Tests and performance 

We tested the set of routines here presented, mainly to confirm the correctness of the returned results. In 
addition, we ran a "performance test" of symgen determining the infinitesimals for Kamke's set of 246 
non-linear 2 nd order examples. Also, although the primary goal of odsolve is to complement the standard 
Maple dsolve in solving non- classifiable ODEs, we ran a "comparison of performances" in solving these 
Kamke examples too. The idea underlying both performance tests was to see the ability of the symgen 
command in determining infinitesimals for a well known set of equations, and to compare the possible 
efficacy of a SM-based solver with that of a classification-based solver such as dsolve. 

5.1 Test of symgen 

As explained in sec.0, the SM approach for solving ODEs involves two main steps: the determination of 
pairs of infinitesimals and their subsequent use in the integration process. The first performance test, thus, 
concerned the explicit determination, by the symgen command, of infinitesimals for the above mentioned 
246 non-linear examples of Kamke's book. To perform these tests we prepared 3 files containing the input 
of all these ODE examples. The table below displays the total number of successes, the average time 
spent with each solved/fail case, and the number of successes of each of symgen's six algorithms when 
the whole set of ODEs was tackled using only one of themfl. 



File 


ODEs 


Symmetries 


Average time 


infinitesimals determined by 


only 1 


many 


ok 


fail 


abacoi 


2 


3 


4 


5 


abaa>2 


1 


100 


49 


24 


7 sec. 


22 sec. 


54 


62 


72 


65 


55 


71 


2 


100 


26 


53 


6 sec. 


22 sec. 


67 


53 


74 


68 


56 


70 


3 


46 


17 


22 


26 sec. 


54 sec. 


29 


20 


30 


27 


23 


31 


Total: 


246 


92 


99 


ss 10 sec. 


« 29 sec. 


150 


138 


176 


160 


134 


171 



Table 1. Kamke's non-linear 2 nd order ODEs: Infinitesimals determined by symgen. 



We here distinguished the cases in which only one symmetry has been found, since in those cases 
the integration procedure has greater chances of only obtaining a reduction of order. The Kamke ODE 
numbers related to the table above are: 



symgen found 


Kamke's book ODE-number 


only 1 symmetry for 92 
ODEs 


11, 12, 15, 17, 20, 21, 22, 23, 24, 25, 26, 28, 31, 32, 40, 45, 46, 47, 48, 49, 54, 58, 64, 65, 66, 
68, 69, 70, 72, 73, 74 , 75, 76, 77, 79, 80, 82, 83, 86, 87, 89, 90, 91, 92, 94, 96, 97, 98, 100, 102, 
103, 105, 106, 118, 119, 120, 121, 122, 129, 132 , 152, 153, 156, 160, 165, 170, 172, 174, 177, 
183, 186, 187, 196, 197, 200, 201, 202, 203, 204, 206, 208, 213, 219, 223, 224, 225, 226, 230, 
231, 238, 241, 242 


many symmetries for 99 
ODEs 


1, 2, 4, 7, 10, 14, 30, 42, 43, 50, 56, 57, 60, 61, 62, 63, 67, 71, 78, 81, 84, 88, 93, 99, 104, 107, 
108, 109, 110, 111, 113, 117, 124, 125, 126, 127, 128, 130, 133, 134, 135, 136, 137, 138, 140, 
141, 143, 146, 150, 151, 154, 155, 157, 158 , 159, 162, 163, 164, 166, 168, 169, 173, 175, 176, 
178, 179, 180, 181, 182, 184, 185, 188, 189, 190, 191, 192, 193, 205, 209, 210, 214, 218, 220, 
221, 222, 227, 228, 229, 232, 233, 234, 236, 237, 239, 240, 243, 244, 245, 246 



Table 2. Kamke 's ODE numbers for which symgen found symmetries: 77%. 

In summary, symgen fails in 23%, finds one symmetry in 37%, and finds more than one symmetry in 
40% of the examples. As a curiosity, although symgen's routines are also prepared to look for symmetries 
not of polynomial or rational type, it is noticeable that most of Kamke's 2 nd order non-linear examples 
have polynomial symmetries. 

8 The infinite. sim.als found by symgen as well as the answers found by odsolve for Kamke's non-linear 2 nd order ODEs 
are available at: ^lttp: / /dft. if.uerj.br/odetools. htm 
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5.2 Test of odsolve's integration schemes and comparison of performances 

As explained in Sec.[I|, for 2 nd order ODEs the knowledge of one or more symmetries may or may not lead 
to the desired answer. For some cases the routines may only obtain a reduction of order. So, to evaluate 
the solving capabilities of the symmetry scheme here presented we ran a performance test for odsolve 
too, and compared the results with those obtained by using the standard Maple dsolve in solving these 
246 Kamke non-linear examples. The results we found can be summarized as follows: 



symgen found 


dsolve solved 


odsolve solved 


odsolve reduced 


many symmetries for: 


99 ODEs 


57 


98 


1 


1 symmetry for: 


92 ODEs 


18 


36 


52 


symmetry for: 


55 ODEs 











Totals: 


246 ODEs 


75 


134 


53 



Table 3. Comparative performance 

Although Kamke's set of ODEs is just a particular one, some conclusions can be drawn from the table 
above. First of all, the standard technique of classifying the ODE according to whether it matches a 
given pattern (for which we know the solving method) did not perform better in any example. On the 
other hand, the implemented symmetry scheme, which tackles the ODEs without making any distinction 
between them, almost doubled dsolve's performance in obtaining closed form solutions, apart from being 
able to return a reduction of order in more 53 examples. Also, odsolve succeeded in integrating to the 
end 36 examples for which only 1 symmetry was found, turning clear the relevance of implementing 
integration alternatives for these cases. Finally, 98 of 99 examples for which symgen obtained more than 
1 symmetry were solved by odsolve to the end, convincing us that the problem is now mainly related to 
finding these symmetries, not their posterior use in the integration process. 

Concerning the cases in which symgen found one or more pair of infinitesimals but odsolve fails in 
integrating or reducing the order of the ODE, the problem was mainly related to: it was not possible 
to find the inverse of the transformation from canonical variables (r, s(r)) to (x,y(x)), or this inverse 
involved RootOfs and integrals. This was the case in Kamke's ODEs 203, 204, 206, 238 and 246. As an 
example, consider Kamke's ODE number 203, 



ode:= ay (y-l)pt-(a-l)(2y-l)(^ 
ax z \ ax 



for which symgen found one symmetry 
> X := symgen ( ode ) ; 



f{x) y (y - 1) ^f-y = 

ax 



(50) 



X 



£ = 0, rj = (y 2 -y) hypergeom ( [1, -], [^—], 

\ a a 



(51) 



The corresponding canonical variables are given by: 
> canoni(X,y(x) ,s(r)) ; 



x, (52) 



ry{x) I 2 a + 1 
s(r) = I (_z 2 - _z) hypergeom([l, -], [ ], 




To reduce the order using these canonical variables the system solved them w.r.t {x,y(x)}, arriving at 
expressions involving the RootOf an integral; even when the change of variables using those expressions 
was successful, due to the complicated resulting algebraic structure, it was impossible to verify the 
reduction of order, and the subroutine gave up with no result. 
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5.3 Test of symgen with the 1 st order Kamke examples 

Although the main purpose of this paper is to present the extension of the package's commands to work 
with 2 nd order ODEs, some improvements have been done in the abilities for solving 1 st order ODEs too. 
An update of the results obtained by symgen for the 1 st order Kamke's examples is as follows: 



File 


ODEs 


Successes 


Average time 


£ and r) can be determined by 


ok 


fail 


abacoi 


2 


3 


4 


5 


6 


abaco2 


1 


90 


63 


14.0 sec. 


7.7 sec. 


33 


14 


21 


25 


32 


24 


35 


2 


100 


88 


2.5 sec. 


1.6 sec. 


27 


37 


48 


47 


40 


45 


44 


3 


99 


86 


1.2 sec. 


8.2 sec. 


34 


73 


56 


54 


50 


1 


38 


4 


89 


80 


3.5 sec. 


16.3 sec. 


41 


40 


47 


47 


56 





44 


5 


89 


77 


3.3 sec. 


15.7 sec. 


14 


7 


68 


72 


69 





38 


Total: 


467 


394 


» 5 sec. 


as 10 sec. 


149 


171 


240 


245 


247 


70 


199 



Table 4. Kamke's (non- quadrature) ODEs, infinitesimals determined by symgen; 86% 



The previous score (see Table 1. in pj) was 75%, and the improvement is mainly due to adjustments in 
'symgen/2', which is now obtaining the infinitesimals for many rational ODEs, and to the introduction 
of a new subroutine, 'symgen/6', for tackling ODEs of Riccati type. The number and class of Kamke's 
1 st order ODEs for which symgen fails in determining a pair of infinitesimals are now given by: 



Class 


Kamke's numbering (1 st order ODE examples) 


rational 


203, 205, 234, 257, 432, 480, 482 


d'Alembert 


388, 390, 405, 406, 409, 410, 430, 479 


Riccati 


22, 25, 110 


Abel 


36, 37, 38, 40, 42, 43, 45, 46, 47, 48, 49, 50, 51, 111, 145, 146, 147, 151, 169, 185, 219, 237, 250, 253, 
265, 269 


NONE 


53, 54, 55, 56, 74, 79, 80, 81, 82, 83, 87, 121, 128, 197, 202, 206, 331, 340, 345, 351, 367, 370, 395, 456, 
460, 461, 493 



Table 5. Kamke 's 1 st order ODEs for which symgen fails: 14% 



6 Computer algebra implementations of Lie symmetry methods 

There are various computer algebra packages for performing Lie symmetry analysis of differential equa- 
tions H [ljj . It is then appropriate to situate ODEtools in the framework of the packages already existing 
in Maple and other computer algebra systems (CAS)0. 

The range of applicability of ODEtools 

To start with, ODEtools is a package for tackling only first and second order ODEs, while some of the 
other symmetry packages, as for example CRACK and PDELIE jl2|, also handle systems of ODEs and 
PDEs of, in principle, arbitrary order. The project is to extend ODEtools' commands to handle high 
order ODEs soon, and concerning PDEs, ODEtools is presently being developed as an ODE-package 
project, but it is fully integrated with the PDEtools package Q . 

Dynamical symmetries and the solving scheme: input=ODE, output=closed- 
form-solution 

Concerning first and second order ODEs, the goals in ODEtools are in some sense a bit different from 
those of other symmetry packages. In ODEtools the focus is mainly in returning closed form solutions 
to ODEs instead of in performing a complete symmetry analysis. Thus, the idea implemented is to look 
for enough symmetries of point or dynamical type as to be able to integrate the ODE - instead of looking 
for all possible point symmetries. Concretely, the difference is that most of the other symmetry packages 
restrict the search to point symmetries^- and are not prepared to build a closed form solution to the 
ODE from the knowledge of these symmetries (see |3|, [l^]). It is worth mentioning here that the search 
for point symmetries is more systematic due to the existence of an overdetermined system of PDEs as 
departure point, but most of second order ODEs don't have point symmetries. On the other hand all 

9 The information present in this section was gathered from j^, |llj, |l2| |l^, |l4| 

10 Exception made by CRACK; restricted capabilities to determine aynamical symmetries are also found in some other 
packages |3l. 



16 



ODEs have infinite symmetries of dynamical type. Also, though the determination of the dimension of 
the symmetry group (done by other packages) before searching for the symmetries may be useful, this is 
of no concrete help when searching for dynamical symmetries since the dimension of the group is infinite. 

Related to these points, the ODEtools command for determining the symmetries (symgen) does not 
restrict the type of the symmetries it looks for, and the ODE-solver (odsolve) has built-in routines for 



using both types of symmetries to directly return a closed form solution for the problem (see sec. 2.2), thus 
allowing the scheme input=ODE, output=answer in a single step as shown in the examples throughout 
this paper. Furthermore, the program is able to proceed with the integration process of a second order 
ODE even when only one symmetry (of any type) was found, by using it to reduce the ODE's order 
and restarting the solving process with this reduced ODE as input. For the case in which the reduced 
ODE cannot be solved, an appropriate scheme for returning reductions of order was implemented. As an 
example of all this, consider 

^--.( d l) 2 -^)(¥]y-oo S ( X )y^ (53) 



dx 2 y \ dx J \dx 

This ODE does not have point symmetries JT^j , and the routines here presented only succeed in finding 
a single dynamical symmetry: 

> symgen (ode ,way=5) ; 

£ _ _ sm(x)y 2 y! + cos(x)y 3 - y 2 - yx 2 ^ 

y 

The differential invariants corresponding to the infinitesimals above are given by: 

[Jo = as, h = sm{x)y + ^] (54) 

y 

and their introduction as variables lead to the following first integral 

1 ( 'hi ; 2 



C 1 = —[-f+sin(x)y') (55) 
y \dx 



with known symmetry (see sec. 4.2). Finding the solution of the above equation is then straightforward. 



In ODEtools, all these steps are run by giving a single input: 
> odsolve (ode) ; 



V 



l + Ci 2 



■ cos(x) + sin(x)Ci + e- c ^ X C 2 + e- Cl X C 2 Ci 



As a second example, let's consider the ODE shown in |l3| illustrating the use of CRACK in solving 
ODEs using symmetry analysis: 

od e: = 3 p 2 h^-5p 2 ( + (5ph-20ph 3 ) ^-20h A + l6h 6 + 4h 2 = (56) 
dp z \dp J dp 

This ODE resulted from an attempt to generalize Weyl's class of solutions of Einstein field equations 
| p5[ . According to |l3|, to solve this ODE in the framework of CRACK and related set of packages 
(LIEPDE, APPLYSYM, etc.) one should: 

1. call LIEPDE to determine the ODE symmetries; 

2. call APPLYSYM with the ODE and its symmetries and run an interactive process to obtain a new 
ODE of reduced order; 

3. use the second symmetry to integrate by hand the reduced ODE of the previous step, to arrive at 
a parametric solution. 
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In the framework of ODEtools, all these steps are compacted into a single call to odsolvep]. Also, many 
integration strategies using the same symmetries in different manners are optionally available. Using 



differential invariants for instance (see sec. 4.1), odsolve directly arrives at an explicit solution for 



(57) 



Eq.([56|) in terms of the roots of a quartic polynomial: 
> odsolve (ode , dif ) ; 

d RootOf (3 _Z*p 2 - Ci 2 (1 + 2 p 2 C 2 ) _Z 2 + 9 p 2 )p 

== ^ - - 

Ci 2 (1 + 2 p 2 C 2 ) (RootOf (3 _Z 4 p 2 - Ci 2 (1 + 2 p 2 C 2 ) -Z 2 + 9 p 2 )) 2 - 18 p 2 

These roots can be made explicit as well, resulting in an explicit answer without RootOf. In this 
example, odsolve spent 5.7 sec calling symgen to determine the symmetries, and 7 sec more to find the 
differential invariants, change variables to reduce Eq.(|5^) to a first order ODE with known symmetries, 
integrate the latter, and return to the original variables. Moreover, in this case the two symmetries 
mentioned in [jl3[ ] are already in normal form (i.e., they form a closed subalgebra), but in other cases, the 
step (j^) mentioned some paragraphs above would not be possible in a direct manner since the reduced 
ODE will not have the same symmetries as the original one. In such a case odsolve will automatically 
either rewrite the symmetries so that they form a closed subalgebra, or determine a third symmetry when 
that is not possible^]. 

These differences do not mean that one package is better than the other - and CRACK is one of the 
most powerful and complete symmetry packages available - but emphasize that the goals and strategy 
of each of the symmetry packages available are not completely equivalent. In CRACK or the Macsyma 
package PDELIE, for instance, there are routines for handling problems completely out of the scope of 
ODEtools as to search for the symmetries of systems of differential equations. 

The symmetries as particular solutions to the determining PDE Eq.(|2|) 

The computational strategy used in ODEtools for finding the symmetries is also different from the 
strategy implemented in most of the other symmetry packages. The difference is that the ODEtools 
routines directly look for particular solutions for the determining PDE Eq. (|2|) , instead of building a set of 
overdetermined PDEs and trying to solve it exactly || . Part of the motivation for this different approach 
is that an overdetermined set of PDEs for the infinitesimals can only be built when one restricts the 
search to point symmetries and this restriction is not present in ODEtools. Another reason is that even 
restricting the search to point symmetries, when a solution exists, this solution usually comes up only 
after completely solving the whole system of PDEs. Of course this requires the appropriate routines for 
"solving systems of linear PDEs" but this is a major problem on its own, and there is as yet no general 
algorithm to integrate such a system Moreover, most of the schemes for tackling systems of PDEs 
will work only if the system of PDEs satisfies a priori some restrictive conditions - e.g., the functional 
dependence on the dependent variables should be of polynomial type or the routines only work with 
rational function coefficients. 

On the other hand, the search for particular solutions to Eq. (|2]) permits searching for point and 
dynamical symmetries in equal footing, by reducing the problem to one of algebraic type or to auxiliary 
first and second order ODEs as explained in our previous work and without requiring routines for 
solving coupled systems of PDEsp|. Moreover, there are no restrictions (e.g., coefficient fields, functional 
dependence etc.) to the type of ODE which can be tackled, and as a matter of fact, at least with Kamke's 
non-linear examples, the average time consumed to solve an ODE is very acceptable (see Sec]|), including 
in that both the determination and the use of the symmetries for building a closed form solution. 

As an example of this, consider the determination of the symmetries for Liouville's ODE 

ode:. ^ +fl(v ) =0 (58) 

ax z \ax I ax 



11 ODEtools includes user-level commands to run each one of these steps as well, see sec. 
12 In these cases, for second order ODEs, a third symmetry can always be determined Jt|. 

13 For the case of point symmetries, and provided that symmetries of polynomial/rational type exist, there is an interesting 
work due to G. Reid and McKinnon for finding particular solutions to the determining system of PDEs by solving auxiliary 

ODEs Fa. 
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where f(x) and g(y) are arbitrary functions - i.e., the ODE has arbitrary dependence in both the indepen- 
dent and dependent variables. Such a problem is out of the scope of CRACK (the ODE is not polynomial 
in the dependent variable), and it is unsolvable using the reduce package SPDE, which doesn't work 
with arbitrary functions ||. It takes 4.8 sec to symgen's subroutines to determine the following two 
particular solutions (in this case two point symmetries) to the determining PDE Eq.(|J): 

> sym := symgen(ode) ; 

[f = 0, 77 = (l + J J 9{v)dy dy^j e -I g( - y)dy ], [£=(l-J e -S f{x)dx dx^j e$ s(x)dx ^ = 0] 
and, in this case, the first of the symmetries above is enough to, in two extra seconds, integrate the 



ODE constructively using canonical coordinates (see sec. 2.2) via 
> odsolve (ode , can); 




-c 2 = o 



First order ODEs 



In the framework of other symmetry packages, Lie methods are usually discarded in the case of first order 
ODEs. The reason is that the determining PDE for the infinitesimals does not split into an overdetermincd 
system of PDEs, and most of the available symmetry packages are based on the splitting of this PDE. 
On the other hand in the framework of ODEtools, the approach presented in Q put Lie methods as a 
concrete alternative for tackling non-classifiable first order ODEs, and showed a performance in solving 
Kamke's first order examples of 85% (see sec. [3J^), representing the highest performance we are aware 
of in solving these examples using symmetries or any other methods. Of course the ODEtools scheme 
will fail with many other examples as well, but until a general solution to the problem is discovered, this 
performance is convincing us that the approach is worthwhile and a computer algebra implementation of 
Lie methods can be a powerful alternative also for first order ODEs. 

Second order linear ODEs 

Another feature of the approach here presented is that it is possible to tackle linear second order ODEs 
using the Lie scheme, while in other symmetry packages these ODEs are excluded from the symmetry 
analysis process (e.g., SPDE), or tackled in a manner which does not represent a concrete advantage. The 
point is that when looking for point symmetries by trying to solve an overdetermined system of PDEs, 
other symmetry packages would look for symmetries of the form [£ = 0,?/ = F(x)], where F(x) is in 
turn a particular solution to the homogeneous ODE. To find this type of symmetry is then so complicate 
as solving the original ODE, and consequently many authors tend to disregard the use of Lie methods 
for linear ODEs. On the other hand, when looking for particular solutions to Eq.(||), these symmetries 
= 0,ry = F(x)] may arise naturally, as in Eq.(||), and symmetries which are not particular solutions of 
the original ODE can be found as well. For example, consider Kamke's second order linear ODE number 
142: 

-)2 



ode := 5 (ax + b) ^ + 8 a^- + c\/ax + by = (59) 
dx z dx 

In this example, a point symmetry which is not of the form [£ = 0, rj — F(x)] is given byf^| : 

14 This symmetry is determined by the internal subroutines as a dynamical symmetry 
[0,F(x,y, Vl )] 

(5ax + 5b)y! + 3 ay 

5=0, V= " 

5 (ax + b) A/5 a 

In the implementation here presented, the additional conversion of dynamical symmetries into point symmetries, when 
possible, is automatic. 
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> sym := symgen(ode) ; 



= (ax + b) 2/5 = 3y 
L * a ' " 5(a S + &) 3/5J 

Now all linear homogeneous ODEs have the symmetry [£ = 0,?7 = y], so, in the case of second order 
ODEs, only one more symmetry as in above is enough to completely solve the problem. Although for 
linear ODEs the standard schemes of rational, exponential, or special function solutions seem to be the 
most appropriate methods, an implementation of Lie's method by looking for particular solutions to 
Eq.(||) appears as an interesting complement since it can handle ODEs with arbitrary functions or solve 
the problem when tricky changes of variables are involved. 

Self-contained environment for first and second order ODEs 

As the last difference mentioned in this summary, the ODEtools package has the intention of being a self 
contained environment for solving and studying first and second order ODEs. For this reason, the project 
includes a compact on-line manual of solving methods for ODEs, implemented here as the odeadvisor 
command, and commands and options playing the role of tools for tackling ODEs which cannot be solved 
using the defaults. 

Among the most relevant of such tools we can mention the dchange program, for performing general 
changes of variables in mathematical expressions; the buildsym program for determining the symmetries 
of an ODE departing from its solution, and the equinv program, for determining the most general ODE 
simultaneously invariant under a given set of point and/or dynamical symmetries. In a typical non- 
directly solvable ODE, one would first try some changes of variables using dchange, and to determine 



the symmetries using the HINT option of symgen as explained in sec. 2.2 15 . When nothing works, 
one can still consider a similar or simpler problem for which the solution is known, try to build the 
corresponding symmetries using the buildsym command, and then try to understand what would be 
the symmetry pattern of the unsolvable ODE. In connection with this, using the equinv command and 
departing from the symmetries returned by buildsym one may be able determine a family of ODEs 
containing the unsolvable ODE. Finally, if the ODE is classifiable, the odeadvisor command may also 
be of help, giving either bibliographical references or some computational hints (e.g., possibly useful 
changes of variables etc.) to study the problem. 



7 Conclusions 

This paper presented a computer algebra implementation of symmetry methods for solving 1 st and 
2 nd order ODEs. This implementation proved to be a valuable tool for tackling ODEs as shown in sec|| 
resulting in the extension to Maple's dsolve we were looking for. Actually, although the implementation 
of the method involves much more complicate operations than the traditional matching-pattern schemes, 
the generality of the families of ODEs which can be concretely solved with it is rather impressive^. 
Moreover, the solver can be used in an interactive manner, that is, one can participate in the solving 



process (the HINT and way=xxx options, the gun variables and the integration schemes; see sec. 2.2), 
achieving in this manner a significant extension of the solving capabilities of the scheme. This possible 
user participation plays a fundamental role, because the functional form for the infinitesimals may be 
out of the scope of symgen's defaults, and because odsolve's defaults, even for the integration schemes 
(see sec. f|), may not be the best choice for all cases. 

Also, two routines, equinv and buildsym, were designed to tackle the inverse problems, that is to 
find the most general ODE simultaneously invariant under a given set of symmetries, and to find the 
symmetries of an unknown ODE when its solution is given. These commands, together with the HINT 
option of odsolve, can be useful in many situations, and for investigating new solving methods. Finally, 
the set of ODEtools commands can be used with educational purposes too, since almost all the relevant 

15 This possibility of participating in the search for the symmetries when desired is a handy and useful feature, almost 
absent in most of the other symmetry packages. 

16 The ODE-solver of ODEtools is now the core of the new ODE-solver of the incomming Maple V Release 5. 
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steps of the symmetry cycle are available as user-level commands, and classification routines for most of 
the 1 st and 2 nd order ODEs are available through the odeadvisor command. 

On the other hand, we also perceived that the complexity of the operations involved in the symme- 
try scheme resulted in a slower solving method if compared with the speed of the standard methods 
based on matching pattern routines. This difference exists when solving 1 st order ODEs (but is almost 
insignificant), becomes more noticeable when tackling 2 nd order ODEs, and turns out to be relevant when 
considering high order ODEs. We then think that an ODE-solver should use both methods, combined 
in an appropriate manner in order to achieve the best performance. Concerning what would be the best 
combination of methods, it is clear that a testing arena should be defined first, and the results will depend 
on this choice. A collection of testing examples most related to real problems seems to be the appropriate 
choice, and in this sense Kamke's book appears to be one of the best candidates available. Then, although 
the natural extension of this work would be to implement symmetry methods for higher order ODEs, we 
think that it is more convenient to invest in fine tuning the merging of symmetry and standard methods 
for the 1 st /2 nd order cases first. This will open the way for more efficient implementations for the high 
order case, as well as permitting us to start more concrete work to tackle PDEs. We are presently working 
in some prototypes in these directions^] and expect to succeed in obtaining reportable results in the near 
future. 
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